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Abstract 

We study the quantum transition from a strongly correlated metal, with 
heavy fermionic quasiparticles, to a metal with commensurate charge or spin 
density wave order. To this end, we introduce and numerically analyze a 
large dimensionality model of Ising spins, in a transverse field, coupled to two 
species of fermions; the analysis borrows heavily from recent progress in the 
solution of the Hubbard model in large dimensions. At low energies, the Ising 
order parameter fluctuations are characterized by the critical exponent zv = 1, 
while above an energy scale, T, there is a crossover to zv = 1/2 criticality. We 
show that r is of the order of the width of the heavy quasiparticle band, and 
can be made arbitrarily small for a correlated metal close to a Mott-Hubbard 
insulator. Therefore, such a correlated metal has a significant intermediate 
energy range of zv = 1 /2 behavior, a single particle spectrum with a narrow 
quasiparticle band, and well-formed analogs of the lower and upper Hubbard 
bands; we suggest that these features are intimately related in general. 
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I. INTRODUCTION 



It has become increasingly clear that studies of magnetic ordering quantum transitions 
in metallic, fermionic systems will be of significant utility in understanding the properties of 
strongly correlated systems like the cuprates or the heavy-fermion compounds — P~3 1 - The 



reader is referred to a recent comprehensive study by Barzykin and Pines |14| comparing 
such ideas with experimental data on the cuprates. 

We begin by discussing some background before introducing the question addressed in 
this paper. Consider a metal in the vicinity of a transition to a ground state with spin or 
charge density wave order. We denote by ^(x, r) the order parameter field as a function of 
space (x) and Matsubara time (r) — the index fi — 1 . . . 3 for the case of vector spin density 
wave order and fi = 1 for the case of scalar charge density wave order. The long distance, 
long-time, zero temperature (T) effective action for obtained after integrating out all the 
fermionic degrees of freedom takes the form P|.|5]P|- |iT| , |r3| 

S = J ' d d qdu\^^ q) u)\ 2 (r + q 2 + ^\uj\ + UJ 2 / c 2 ) + .... (1.1) 

We have Fourier transformed to wavevector q and Matsubara frequency u, and introduced 
phenomenological constants r, 7 and Cg; the ellipsis represent higher order couplings between 
the ^ fl . We are considering here only the case originally considered by Hertz |I[ in which || 
the ordering wavevector, Q, is not a spanning wavevector of the Fermi surface, and it is 
possible to connect at least two points on the Fermi surface by Q — for a spherical Fermi 
surface this last condition is \Q\ < 2kp, where kp is the Fermi wavevector. In the language 
of a recent paper by one of us, Chubukov and Sokol [13| , we are considering only transitions 
of type B. Under these conditions, the density of particle-hole excitations with wavevector 
Q are linear in energy at small energies, and this is responsible for the dissipative \uo\ term in 
S. The q 2 and u 2 p|-[Tl|,p^1 terms are simply regular terms controlling spatial and dynamic 
fluctuations. Alternatively, the u 2 term may be viewed as the simplest addition that makes 
the order parameter correlator have the correct 1/uo 2 decay at large frequencies. 

In mean-field theory, 5* undergoes an ordering transition when the Landau-like parameter 
r = — hence r > measures the distance of the metal from the ordered state (we caution 
that an action like 5* is not appropriate for r < Jl3|). For the usual reasons, the mean- 
field transition has a correlation length exponent v = 1/2. At low energies, the \u\ is the 
most important term controlling the energy scale — comparing this with the q 2 term we get 
a mean-field value for the dynamic critical exponent of z = 2 M . Later we will consider the 
solution of a model for such a transition in the limit of large dimensionality, where z and v 
are not separately defined, but the energy scale exponent zv is. Hence we prefer to quote 
the value of zv which is zv — 1 for the present mean-field theory. It has been argued on 
theoretical and phenomenological grounds JT^J^J^-|n],|TB[ that, under suitable conditions, the 
nearly-ordered metallic state should exhibit a crossover at higher frequency or temperature 
scales to region with the effective exponent z — 1, or in mean-field theory, zv — 1/2 (at the 
z — 1 fixed point we still have v — 1/2 in large d mean-field theory, but v w 0.7 in d — 2). 
It has also been suggested that this z = 1 behavior is observed in the cuprates not only at 
very low doping H but also at moderate dopings |6|,|T4f. In a simple mean-field analysis of 
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S, the energy/temperature scale above which zv = 1/2 behavior appears is of order 7Cq, as 
that is when the uo 2 term starts becoming more important than the |u;| term. 

The question we address in this paper is the following: under what conditions does this 
crossover from zv = 1 to zv = 1/2 actually occur ? In other words, is there a reasonable 
scenario in which the crossover scale ~ jCq is significantly smaller than all other higher 
energy cutoffs of the critical behavior, so that zv — 1/2 behavior is clearly observable ? At 
sufficiently high energies the behavior of the system must become dominated by lattice scale 
cutoffs, and hence non- universal, and it is therefore important to have the cutoff ~ 7Cq be 
lower than such cutoffs. 

We will answer this question here using the solution of a large dimensionality model of 
spinless fermions which has a charge density wave-like order parameter with a Z 2 symmetry 
(the solution uses recent progress in understanding the Hubbard model in large dimen- 
sions [16-20fl). Our model exhibits a transition between metallic states with and without 
a mean value for the Z 2 order parameter. We will show that there is indeed a reasonable 
scenario in which the vicinity of this transition exhibits zv = 1/2 behavior over a signifi- 
cant intermediate energy scale. Although we will explicitly display results only for a specific 
model, our arguments are quite generic; we expect that similar results apply to other models 
of ordering transitions in fermionic systems in infinite dimensions, including those with order 
parameters with a larger symmetry e.g. the 0(3) symmetry in spin density wave transitions. 

We now turn to a qualitative discussion of the region of parameter space where such 
behavior occurs. Rather than introducing the infinite-dimensional model at this point, let 
us discuss a more familiar example in two dimensions; this example has the virtue of having 
phases closely analogous to all of the phases of the infinite-dimensional model (and perhaps 
more). This will make the physical meaning of the phases transparent, allowing us, then, to 
jump directly to the phase diagram of the infinite-dimensional model. 

Consider spin- 1/2 electrons moving on a square lattice with first (ti) and second (t 2 ) 
neighbor hopping matrix elements and short range repulsive interactions of scale U which 
prefer spin density wave ordering at the wavevector Q = (71", ir) (See e.g. [jlvj). We need 
second neighbor hopping to avoid effects from nested Fermi surfaces. On general grounds 
we anticipate 4 distinct ground states for this model at half-filling: 

(A) Metal — this is a Fermi liquid whose Fermi surface crosses the boundary of the magnetic 
Brillouin zone. Such a Fermi surface is not nested and is also of the type in which spin 
fluctuations are described by an action like S. 

(B) Metal with spin density wave order — the magnetic order causes gaps over portions of 
the Fermi surface, but the system remains metallic. 

Some details of the evolution of the Fermi surface and critical properties of the transition 



between phases (A) and (B) were discussed in Ref [[yj. At very large values of U we can 
have two additional insulating phases 

(C) Insulator with Neel order — all charged excitations are now gapped. The appropriate 
model for spin excitations is the spin-1/2 Heisenberg spin model with first and second 
neighbor exchange — the so-called J1-J2 model. 

(D) Insulating quantum paramagnet — or a Mott-Hubbard insulator. The J1-J2 model can 



have a ground state in which the magnetic order is short-ranged 0,0] • This state can also 
have spin Peierls order f2~2| , but this is incidental to our considerations here. 
Determining, in this two-dimensional model, the topology of the phase diagram with phases 
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A, B, C, D, and possibly others, the nature of the phase transitions, and of the possible 
multicritical points, is a formidable problem which we shall not address here. 

Turning then to the large dimensionality model to be introduced later in this paper, we 
show in Fig [l] its phase diagram. The phases are closely analogous to phases A, B, C, and 
D above and are labeled as such in the figure. The model has a charge density wave-like Z 2 
order parameter (r z ), but its role is analogous to the Neel order parameter above. 

We are now ready to state the main result of this paper. We are interested in critical 
properties of the transition from phase A to phase B. Consider the region of phase A which 
is close to the Mott-Hubbard insulator D. We identify this region as a 'strongly correlated 
metal', and has been shaded as such in the Fig [I]. In such a metal p!9|j20|j23|j2^] the fermionic 
quasiparticles are quite heavy and form a narrow band of width T around the chemical 
potential. A large part of the single particle spectral weight is in the higher energy lower 
and upper Hubbard bands. Now consider the ordering transition from the strongly correlated 
metal to the metallic, ordered, phase B; in other words, the transition from A to B, not too 
far from the multicritical point M. We shall show that in the large dimensionality limit, the 
order parameter fluctuations are characterized by zv = 1 only at energy scales below T, 
while zv = 1/2 behavior takes over at the intermediate energy scales between V and y/TU . 
For all energies below y/TU, and in finite, but large, dimensions, the action S in Eqn. ( |1 . 1| ) 
is an excellent starting point for describing order parameter fluctuations and the crossover 
between zv = 1 and zv = 1/2 criticality, as examined in Ref [ 13| . As T can be made small 
by moving towards phase D, there is a significant window of intermediate energies with 
zv = 1/2 order parameter fluctuations. (For extremely small values of T, it may also be 
necessary to consider the energy range between y/TU and U — this regime also has zv = 1/2 
criticality, but with a significant non-universal renormalization of the parameter cq in S 
at the scale y/TU). These results suggest an intimate relationship between a low energy 
crossover to zv = 1/2 behavior, heavy quasiparticle bands, and the removal of spectral 
weight to the lower and upper Hubbard bands — we conjecture that this relationship is more 
generally valid. 

Finally, we mention that a recent study [^5| of liquid 3 He considers a large dimensionality 
model whose phase diagram has many similarities to Fig [I], and which also appeals to the 
proximity to a multicritical point closely anologous to M. 

The infinite dimensional model will be introduced in Section [TI| and its ground states 
discussed in Section |H|. The main conclusions will be reiterated in Section [TV. 



II. THE MODEL 

We motivate our model from the d = 2 spin-1/2 model considered in Section |. In the 
vicinity of the transition between phases A and B, we need only focus on the spin density 
wave order parameter and portions of the fermi surface which are close to pairs of points 
that can be connected by the ordering wave-vector Q. Let k\ and k 2 = k\ + Q be one such 
pair of points; in general, there will be other pairs of points, usually related to k%, k 2 by 
symmetry operations of the square lattice — these can be treated in a similar manner and 
are not considered explicitly. We introduce two species of fermions ci a , c 2a (a =T, I is the 
spin index) representing Fourier components of the electron in the vicinity of k\ and k 2 . The 
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low energy effective action for the vicinity of the boundary between A and B can be written 
as [Il3l 



S = J d d xdr [(V^) 2 + (d^rf/cl + f^l + c\ a (d T - v, ■ V)c lQ + c\ a (d T - v 2 ■ V)c 2o 

-A*,(cL<^ + H.c.)], (2.1) 

where vi, v 2 are the Fermi velocities at fci, k 2 , and <t m are the Pauli matrices. After integrating 
out the fermions, the effective action for VP^ should be of the form ( |1 . 1|) . The dissipative 
term is of the form as long as V\ and v 2 are not anti-parallel i.e. there is no nesting. 

We now abstract from S the essential ingredients for a simple model of the ordering 
transition in large dimensions. First, we reduce the symmetry of the order parameter from 
0(3) down to Z 2 (the symmetry of a charge density wave order parameter) by discarding 
the spin index a. Then, we replace the scalar field \I/ by an Ising spin t z — its 'kinetic' term 
then becomes a transverse field, g. (For the model with 0(3) symmetry, the Ising spins in 
a tranverse field would be replaced by 0(3) quantum rotors). Finally, we place the system 
on a lattice (not the same lattice for which S was the continuum limit) and write down the 
Hamiltonian 

H = -9 ^2 Txi ~~Z ^ ^ TziTz j ~ ^ XX C li Cl * + C 2i C 2i) -7= ^(hijCiiCij + t 2ij C 2i C 2j ) 

i>j i * ij 

* T zi( C li C 2i + c ] 2i cu). (2.2) 



We have introduced the two Pauli matrices r z , t x , an Ising spin exchange J^, a chemical 
potential /z, hopping matrix elements tuj and t 2 ij and a fermion-Ising coupling A. The 
parameter Z has different meanings depending upon the particular large dimensionality 
limit chosen Hl6,[l9,p0| — for a regular d dimensional lattice, Z = d; for the Bethe lattice, Z 



is the co-ordination number; and for a fully connected cluster with random hopping tuj, hij 
between all pairs of site, Z is the total number of sites. The Hamiltonian H contains all low 
order terms consistent with the Z 2 symmetry 

T z -> T X T Z T X = -T z ; T x -> T X T X T X = T x ] C X C t ] C 2 -> ~C 2 . (2.3) 

In the context of the original lattice model of electrons, this Z 2 symmetry is simply the 
sublattice interchange symmetry which is broken by a charge density wave ordered state. 

We will study in this paper the Z = oo limit of H. Further, we will focus only on solutions 
of the mean field equations in which all sites are equivalent. This amounts to neglecting 
possibilities in which there is a staggered ordering of the Ising spins r z . Such solutions do 
occur, and are often lower in energy than the ones we consider [^,^,^|; however we are 
not interested in them on physical grounds. One can also appeal to the fully connected 
random clusters for which the spatial uniform solutions are expected to be true ground 
states |3|lg|23|j2§. 



We will make an additional, final, simplification. We will work with models in which 
tuj = t 2 ij. Although there is no symmetry enforcing this equality, in finite-dimensional, reg- 
ular lattice models such an assumption would drastically change the physics of the transition. 
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This is because, as will become clear below, for tuj = t 2 ij there are additional conserved 
charges which have a significant effect on the critical theory [[j]. However at Z = oo, all 
of the non-trivial critical behavior is in local correlators which are expected to behave in a 
similar manner for tuj = t 2i j and tuj ^ t 2i j. The numerical analysis required to solve the 
mean-field equations becomes much simpler at tuj = t 2 ij- 

We conclude this section by writing down the mean-field equations. The equations take 
their simplest form after performing a rotation of the fermionic fields 

c a = ^(c 1 + c 2 ) C b = -j=(d -ca). (2.4) 



Under the Z 2 symmetry ( |2.3Q we have c a — > Cb and c& — > c a ; this makes it clear that we 
can simply think of the a and b fermions as moving predominantly on the two sublattices 
of the model underlying S. For Z = oo, the model is mapped onto a single-site mean-field 
hamiltonian supplemented by a self-consistency condition. Following Refs [[17],|19|,^(|] , for 
t\ij = hij = t for nearest neighbors on the Bethe lattice, or for the fully connected cluster 
with tuj = t 2 ij random but Jjj uniform, the single site mean-field Hamiltonian is 

H M f = ~gr x - M4 c a + c\c b ) - Jmr z - Xr z {c\c a - c\c b ) 

+ (takfLfak + ebkftjbk ~ V ak c\f ak - H.c. - Vbkclfbk ~ H.c.) . (2.5) 

k 

We have introduced fermions /„, /& to model the environment of the site of interest. The 
couplings of these fermions and the parameter m are determined by the self-consistency 



conditions 17.20,19.23 





m 


= (r.) 


t 2 G a 




_ \- \V a k 2 
^ iu - tak 


t 2 G b 




\Vbk 2 



(2.6) 

where G a is the Fourier transform of the a fermion Greens function — ^c a (0)<4(r)^) and 
likewise for Gb- 

Note that the total number of a fermions and b fermions is separately conserved. This is a 
consequence of the choice tuj = t 2 ij] we reiterate that while such a choice and the additional 
conserved quantities would be dangerous in the finite dimensional theory, it's effect on local 
correlators in the Z = oo limit is expected to be innocuous. 



III. GROUND STATES OF THE MODEL 

This section will consider solution of the Z = oo model Hmf in Eqn (|2.5| ) along with 
the self consistency conditions ( |2.6| ). We will consider subsection [111 A| the exact solution in 



the atomic limit t = 0, J = 0, followed by a numerical solution of the t ^ 0, J = case in 
subsection [111 Bl The most general t ^= 0, J > case will be considered in subsection [ill C 
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A. Atomic Limit, t = 0, J = 



All of the sites are now independent. Each site can have either 0, 1, or 2 fermions and 
the eigenenergies of H can be easily determined. They are 

Eq = ±g , 2 states 

Ei = — p ± \J X 2 + g 2 , both doubly degenerate, 4 states 

E 2 = -2p±g , 2 states. (3.1) 

A key observation is that if we take the lower eigenvalue for each particle number, the states 
map exactly onto those of the Hubbard model with U/2 = y/X 2 + g 2 — g. Moving away from 
the atomic limit, the mapping of the low energy states to the Hubbard model will continue to 
hold as long as t, J < g. This mapping will be very useful to us in the subsequent discussion. 



B. J = 

We will study the ground states at fixed values of A and t as a function of g. This 
corresponds to the y-axis in Fig [l] which presents results at t = 1 and A = 2.5. 
First consider the limits of large and small g. 

For large g, the Ising spin flips rapidly between its up and down states as a consequence 
of the gr x term. The value of r 2 averages out to zero, and the fermions effectively do not 
see the Ising spin. The fermion spectral function p a ,b(ty — ImGa^fi) is then simply the 
semi-circular density of states of free fermions: 



p a (Q) = Pb {Q) = — ^At 2 -(n-/2) 2 ; g^,oo (3.2) 

for |fl — p\ < 2t and zero otherwise. 

The behavior at small g is a little more subtle. The Ising spin now fluctuates slowly 
between its up and down states and the fermions have plenty of time to respond to its 
instantaneous orientation: this yields fermion bands centered around ±A. As the Ising spin 
is equally likely to be up or down, the fermion spectral function is an equal superposition of 
the two possibilities 



Pa (Q) = Pb (tt) « — 



At 2 - (fi - A - p) 2 + J At 2 -(n + x-p) 



g small (3.3) 



where it is assumed that the square roots vanish when their arguments are negative. For 
A > 2t there is a window of energies where ( |3.3D predicts a gap; in reality, we will only 
have a pseudogap whose origin is related to the phenomenology discussed by Kampf and 
Schrieffer ||28|| . For g > t (see Section [III A|) it is more appropriate to think of the pseudogap 



as arising from the formation of the 'upper and lower Hubbard bands'; there is, however, 
no fundamental distinction between the two mechanisms and they continuously evolve into 
each other. 

One might also wonder why, for small g, we do not expect the Z 2 symmetry to be 
broken, thus causing the Ising spin to pick a definite orientation and p a ^ Pb- This, however, 
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does not happen for the same reason that the insulating solution of the Hubbard model in 
Refs [19|20],^3| is not ferromagnetic. One can test instability towards 'ferromagnetism' by 
measuring the response to an external field — even though the local susceptibility is infinite in 
the insulating solution, there is an internal field from the neighboring fermions which causes 



the net local effective field to vanish [23], and the uniform susceptibility remains finite; the 



reader is referred to the discussion in Refs |^|^] for more discussion on this important 
point in the context of the Hubbard model. 

We have numerically studied the crossover from small to large g by solving the Eqns ( |2.6| ) 
using the exact diagonalization method |3"0"| , |3"T| ] . We used as many as 8 sites each for the a 
and b fermions and worked at half-filling. The results were very similar to those obtained 
in the Hubbard model |19| , ^0| , |23[| . We obtained two classes of solutions, a metallic solution 
present for g > g c2 and an insulating solution present for g < g c \ with g c \ > g c2 . For A = 2.5 
and t = 1 (the parameters used in Fig [I]), the energies of the two solutions crossed each 
other at a value of g ~ 0.7 which was quite close to g = g c2 . We are therefore uncertain 
as to whether the metal-insulator transition was first or second order. However this issue 
is really peripheral to what we are interested in here, and so we did not do the necessary 
investigation of numerical fine structure to resolve it. 

What we do care about is the behavior of the metallic solution for g close to but greater 
than g c2 — i. e the shaded region of the strongly correlated metal in Fig. |l|. In this region, it 
is appropriate to use the continuous transition at g = g c2 to obtain a scaling description of 
the response functions. There is a low energy scale, which we denote by T, which vanishes 



as g approaches g c2 as p3 



9 - 9c2- 



(3.4) 



The scale T controls the width of the quasiparticle band (~ T) or the effective mass (~ 1/r). 
Our main interest here is in the behavior of dynamic correlations of the order parameter r z . 
Its scaling properties can be deduced from the elegant critical theory of the Mott transition 
provided recently by Moeller et. al. [E3J. (Our model H has only a Z 2 symmetry and so the 
critical theory will be a self-consistent Kondo-like model |24| but with a planar anisotropy 
in the exchange constants). For the response, xioc, to a local field coupling to t z 



Xloc{UJ) 



dre*" (t z (0)t 2 (t)) , 



we have [241: 



X 2 u? + ... 



M < A, 



(3.5) 



(3.6) 



with cf) a universal scaling function, A is an upper cutoff of order A or g, and constants X±, X 2 
contain the corrections to scaling contributions of the higher energy excitations associated 
with lower and upper Hubbard bands. The latter excitations have an energy ~ A and 
hence their contribution can be expanded in a smooth power series in even powers of u. On 
dimensional grounds we expect X\ ~ 1/A and X 2 ~ 1/A 3 . 

We show in Fig |2| a test of the scaling form ( |3.6| ) using the numerically computed value 
of xioc at four values of g > g c2 — the collapse of the data onto a single scaling curve is quite 
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reasonable. The scaling function can be chosen such that 0(0) = 1 and must satisfy the 
asymptotic limits 

/ / n = f 1 - ci|x| for |x| -> , s 

1 c 2 /:r 2 for |x| — > oo ' 

for some positive constants c.\,c 2 . The non-analytic |x| behavior at small frequencies is 
a consequence of the damping of spin excitations from the finite density of states at the 
Fermi level in the metal, while the x 2 behavior at large frequencies follows simply from the 
spectral representation of the response functions in the critical theory. These asymptotic 
limits suggest the simple interpolation form 

1 + C\ \X \ + C 2 X Z 



which satisfies ( p.7| ). We also show in Fig. ^| a fit of the measured scaling functions to 
It is apparent that the interpolation form works rather well, and that the accuracy of the 
numerical results is not sufficient to distinguish between the true and the approximate 



form (3.8). For the parameters chosen, the best fit values were c\ ~ 0.07 and c 2 = 0.32; the 
value of c 2 is quite reliable, but the same cannot be said of c\ — slightly different choices 
in the fitting process gave values of c\ differing by a factor of 2, but only a few percent 
changes in c 2 . However, C\ always remained significantly smaller than c 2 — indeed, the linear 

dependence of is barely visible in Fig|2|. We had to compute xioc at fairly large values of 
g (g > 2.0), lying well within phase A, before the metallic |u;| behavior was clearly evident. 

The scaling behavior of the q = susceptibility, Xq=o(w), which determines the dynamic 
response of the t z to a spatially- uniform but time- dependent field acting on the r z , is some- 
what more complicated (Note that since t z mimics the density wave order parameter at 
a wavevector Q in the physical system, the q = susceptibility of H is really the phys- 
ical staggered susceptibility). Unlike, xioc this susceptibility is finite at the transition at 
g = g c2 [|T9|j20| , |23[1 ; as a result the interesting frequency dependence of Xq=o is really in a 
correction to scaling contribution. With the limited accuracy of the numerical method we 
are using, it would be quite difficult to obtain the scaling results from the data. We there- 
fore restrict ourselves here to a qualitative discussion, which will be sufficient to extract the 
physics we are interested in. Because of the similarity in the low energy structure of H to 
the Hubbard model (in particular, the two-fold degeneracy of the lowest energy state E% in 
( ^.1[ ) in atomic limit), we expect that Xg=o wm behave in a manner similar to the uniform 



spin susceptibility near the Mott transition. From arguments similar to those in Ref |23] we 
may then deduce that Xq=o behaves roughly like 

Xqlo ~ X + caxTol, (3.9) 

where the constant X ~ A and C3 is a dimensionless constant of order unity (on the hyper- 
cubic lattice, the above equation actually requires tuj 7^ faij)- Combining ( |3.6|) , (|3.7p and 



9|), we can now deduce the frequency dependency of the uniform dynamic susceptibility: 

-if \ v , j c 3 (r + ci\u\ + c 2 uj 2 /T) for \u\ < VTA , , 

X - oM ~ X + j 03(1/^ + X 2 u 2 /Xl + ...) for VFA < M < A ■ (3 - 10) 

The first of these results in the energy range \oj\ < VTA will be used below to reach the 
main conclusions of this paper. 
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C. J > 



The exchange coupling J plays a very simple role in the Z = oo mean-field theory. It 
is apparent from Eqns ( |2.6|) that in phases with m = (t z ) = the single particle Green's 
functions are in fact independent of J. The results at J = thus continue to apply for a 
finite range of values of J. 

There is however a change in the t z correlation functions. It is not difficult to show, 
using an argument quite similar to that in the original Curie- Weiss mean field theory, that 



xiM -J. (3.11) 



j=o 



The finite J uniform susceptibility therefore diverges at J = Xqloi 00 ) 3 Q 5 which is the point 
we have onset of a nonzero (r z ). The phase boundary between phases A and B and between 
D and C in Fig [1] was determined in this manner. The phase boundary between phases C 
and D requires computations in the phase with (t z ) ^ 0; as this phase boundary is of no 
interest to us here, we did not carry out the rather involved computations required — the 
boundary between C and D, shown in Fig |I| is just an educated guess. 

Fig [I] also shows a multicritical point M, which is the point where all four phases would 
meet if the metal-insulator transition was second-order. The correlation functions on the 
( T z) — side of this critical point are however simply related to those at the g = g C 2, J = 
critical point discussed in Section |III B| : the single particle correlators are unchanged, while 
the relationship in the spin susceptibility follows from ( |3.11| ). 

The central interest of this paper is in the nature of the order parameter fluctuations 
near the phase boundary between A and B, at a point not too far from M (Fig [l]). From 
( j3.11| ) and (|3.10 ) we see that the static uniform susceptibility diverges at J = X, and the 



Landau parameter r = X + T — J. The nature of the dynamic susceptibility near the phase 
boundary between A and B also follows from ( |3.11| ) and ( |3.10| ). From these equations we 



see that the infinite dimensional transition has zv — 1 at the lowest energies; however this 
behavior is only present for energies smaller than T, and zv — 1/2 criticality takes over for 
larger energy scales. Observe, however, by comparing the two equations in ( [3.10D , that there 
is a significant change in the coefficient of the u 2 term at a scale VTA. 

We can in fact also make a crude connection between the dynamic susceptibility com- 
puted here and the finite-dimensional action S for discussed in Section |. One perspective on 
the meaning of the results for dynamic correlators in the infinite dimensional model is that 
these specify the input form of the effective action that must then be used to understand 
fluctuations in finite dimensions. With this point of view, our results allow us to deduce the 
following effective action for the order parameter ^ in finite dimensions: 

S = /dVM^(9,^)iV+X g "ioH) + ---> (3-12) 

where the q 2 has been added on phenomenological grounds and, as in (|1 . 1| ) , the ellipsis 
represent non-linearities not explicitly displayed. From ( P-lOj ) and ( |3.11| ) we see again that 
r = X+T—J, and that the two actions (|1.1|) and (|3.12|) have identical frequency dependencies 
in the frequency range |a>| < VTA where 
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S = Jd d qdu\^ il (q 1 u)\ 2 {r + q 2 + c^c 1 \u\ + c^c 2 uJ 2 /Y) + .... (3.13) 

The insight gained from the present approach is that we now have an estimate of the scale 
at which the crossover from zv = 1 to zv = 1/2 behavior occurs — it is the small energy T, 
as had been claimed earlier in the paper. The crossover itself is described by the action S 
in ( |3.13 ), which is an excellent approximation for energies below vTA. In the higher energy 



range VTA < \u>\ < A we should use instead the second result for x q =o m ( |3 . 1 0|) — this range 



still has zv — 1/2 behavior, but there is a non- universal crossover in the coefficient of the u 2 
term at the boundary \oj\ ~ VTA. So if F is extremely small, this second regime of zv — 1/2 
behavior will be the most important, and the very low energy regime of zv = 1 behavior 
will be unobservable. 

Before concluding, we also note that a similar approach can be used to describe the 
ordering transition between the insulating phase D and C. In this case T = and we are 
always in the second of the regimes in (|3.10|) — the transition therefore has zv — 1/2 down 
to the lowest energy scales. 



IV. CONCLUSIONS 

This paper has examined a simple, infinite-dimensional model for spin or charge density 
wave ordering transitions in strongly correlated metals. We were interested in the case, 
originally considered by Hertz [Q , in which Fermi surface geometry was such as to induce a 
|a; | damping in effective action for the critical modes of the order parameter. As a result, 
the order parameter fluctuations have zv = 1 at the lowest energy scales. We examined 
the conditions under which there is a well-defined crossover at intermediate energy scales 
to zv — 1/2 criticality. We showed that, under suitable conditions, it was indeed possible 
to have a low energy scale T such that the zv = 1 behavior was restricted to \u\ < T, and 
universal zv — 1/2 behavior appeared in the energy range r < \lu\ < vTA; here A is an upper 
cutoff of order the repulsion energy between the electrons — A ~ U. For all energies smaller 
than a/TA the order parameter fluctuations in large, but finite, dimensions are expected to 
be well described by the action S in Eqn (|1.1|) or (|3.13|) : a detailed analysis of the finite 
temperature crossovers and fluctuations in a d = 2 system controlled by S has already been 
presented in Ref ||13|| . For frequencies larger than vTA, the order parameter fluctuations 
still have zv = 1/2 — the main caveat to keep in mind is that there is non-universal change in 
the value of Cq (see (|1 . 1|) around uo ~ vTA which is not describable by the simple frequency 
dependence in S. 

The key feature of the above scenario is of course the presence of the low energy scale V 
in a strongly correlated metal. In infinite dimensions, V vanishes at the transition to a Mott 
Hubbard insulator, which we assume (see Fig [1]) is in the vicinity of the strongly correlated 



metal. In this same region [[T£],^(],[^| , the single-particle spectral function has a narrow 
quasiparticle band of width T, and additional spectral weight at the energies ~ A which 
form the analogs of the upper and lower Hubbard bands. This window of energy between 
the quasiparticle excitations and the localized Hubbard excitations is directly responsible 
for the zv = 1/2 criticality in the order parameter fluctuations. 
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FIGURES 



FIG. 1. Phase diagram of the Hamiltonian H ( [2. 2D at half-flung in the Z = oo limit. We 
have chosen t = 1 and A = 2.5. The metal-insulator transition between phase A and D is assumed 
to be a second-order transition at g = g C 2- However, the phase D continues to be metastable for 
9ci < 9 < 9c2 and the energy of this state crosses that of A at a point very close to g = g C 2- 
Thus this transition could be weakly first-order, in which case g C 2 and the multicritical point M 
will be in a metastable region. However, the critical theories at g = g c i and M could still be 
used to describe the shaded strongly correlated metal. The Ising order parameter is (t z ). We are 
interested mainly in the Ising transition between A and B from the shaded region of A: in this case 
the Ising fluctuations have zv = 1 at low energies and zv = 1/2 at intermediate energies. The Ising 
transition between D and C has zv = 1/2 at low and intermediate energies. The phase boundary 
between the B and C is of no real interest to the discussion in this paper, and its position is based 
on an educated guess; positions of other phase boundaries were numerically determined. 

FIG. 2. Scaling plot of the local Ising spin correlation function xioc f° r i = 1, A = 2.5 in phase 
A of Fig p]. The values of T with their respective values of g are: g = 0.75, T = 0.0056; g = 0.8, 
r = 0.0090; g = 0.85, V = 0.0128; g = 0.0, T = 0.0179. The fit function <j) is defined in Eqn flpj), 
and the values of the fitting parameters are quoted below it. 
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